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Fulde-Ferrell-Larkin-Ovchinnikov oscillations in the pair-pair correlations. For 
m = mi a, kinetic energy driven macroscopic phase separation into low-density 
superconducting domains and high-density polarized walls takes place. For m > 
1712 the domains fully localize, and the system eventually becomes a ferrimagnetic 
insulator. 
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1. Introduction 

In the century-long search for a mechanism which unveils the origin of 
superconductivity in so many different materials, a key ingredient is to understand 
the role of magnetic correlations in the formation of the superconducting (SC) pairs. 
Apart form the Meissner effect, different magnetic effects have been described in 
superconductors, such as the ferromagnetic to superconductor transition in heavy 
fermions compounds[lJ, or the magnetic field driven SC-insulator transitions in 
some two-dimensional high- Tc samples [2j. The very recent discovery of high- Tc 
superconductivity in iron based layered pnictides [3j , and the observation of coexisting 
micro- and/or mesoscopic SC and magnetically ordered domains[4, 5j which are 
reminiscent of the stripes characteristic of cuprates, confirms the suggestion that 
superconducting and magnetic order are intimately related, and can be crucial for 
each other's stability [6j. On the theoretical side, it has also been predicted that 
SC correlations in presence of non-zero magnetization can rearrange their spatial 
modulation exhibiting Fulde-Ferrell-Larkin-Ovchinnikov(FFLO) oscillations [7J, which 
are expected in case of strongly anisotropic (for instance one-dimensional) systems. 
It is a remarkable very recent achievement their observation in both heavy fermions 
systems and one dimensional Fermi gases [9, 8j. 

It is generally accepted that the above behaviours are to be ascribed to the 
presence of electronic correlations in the different systems. The reference model to deal 
with electron-electron interaction on lattices is given by the Hubbard Hamiltonian, in 
which the main contribution of the Coulomb repulsion to the model Hamiltonian is 
identified with the on-site interaction of electrons with opposite spins {U term). The 
presence of a SC phase at /7 > for this model is still matter of debate, though 
some encouraging results have been achieved (see for instance [lOj and the discussion 
therein). Many extensions of the Hubbard model have been proposed, based on the 
inclusion of the first contributions to the Coulomb repulsion other than the U term. 
Among them, a generalization motivated independently by Hirsch[llJ and Gammel 
and Campbell [12j amounts to include the nearest neghbors off-diagonal interaction 
X, which turns out to modulate the hopping term depending on the occupations of 
the two sites. In this case, the model Hamiltonian reads 

H = - ^ [l-^(nia+^ja)]4a^i^+^X^^^t^ii+^X](^*t-^^i) (1) 

<ij>a i i 

where creates a fermion with spin a = {t,i}, ^ denoting the opposite of a, 
^icT = cJ^Qcr is the <j-electron charge and (ij) stands for nearest-neighboring sites. The 
parameters u and x are the dominant diagonal and off-diagonal contribution coming 
from Coulomb repulsion, and the lower case symbols denote that these coefficients have 
been normalized in units of the hopping amplitude. Here we have also included an 
external magnetic field h. Moreover, N is the number of electrons on the d-dimensional 
L-sites lattice, so that n = N/L is the average filling value per site. The model has 
been extensively studied in the literature at = [13j. In particular, since H is not 
invariant under particle-hole transform, it has been proposed to model a theory of 
hole superconductivity, [llj 

While the effect of the off-diagonal Coulomb repulsion (also known as bond-charge 
interaction) has long be disregarded since it is usually smaller than the diagonal 
contributions of Coulomb interaction (Hubbard and extended Hubbard models), it 
has by now become clear that what matter is instead its amplitude compared with 
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that of the hopping term. In fact already for x > 1/2 a scenario quite different from 
that of the Hubbard hmit emerges [T^ [T5] even in d = 1. In particular, for not too 
large on-site Coulomb repulsion u < Uc{x) a new metallic phase appears, characterized 
by dominant SC correlations and nanoscale phase separation [16j (NPS). The latter 
amounts to the coexistence of two conducting phases with different densities [17J, and 
is driven by the charge degrees of freedom. The actual size of the coexisting domains 
turns out to be microscopic at = 0, and to increase in relation to the imbalance 
in spin orientation[l8], suggesting an interplay of charge and spin degrees of freedom 
in the ground state. Since none of the two coexisting phases is separately SC, it is 
reasonable to interpret this emergent SC behaviour as a consequence of the microscopic 
domain size induced, in turn, by spin rearrangement; we are in presence of magnetic 
mediated superconductivity. A further confirmation of this hypothesis comes from a 
very recent result [19j which shows that the critical line for the SC transition can be 
recovered exactly (see figure [TJ left panel) if assuming short range antiferromagnetic 
correlations between just the single electrons. It must be noticed that the appearance 
of antiferromagnetic and metallic behaviour at half-filling, in a regime of moderate 
interaction, as well as its interplay with SC properties, is somehow reminiscent of the 
physics characterizing iron pnictides in higher dimension. 

To better understand the role of magnetic correlations in ground-state of 
Hamiltonian ([T]) and their interplay with the charge degrees of freedom, here we 
shall investigate -through a detailed density-matrix renormalization group (DMRG) 
numerical study [20j- the consequences of the inclusion of the magnetic field as to 
the onset of the SC phase. In principle, increasing the magnetic field should force the 
spin degrees of freedom to align to the external field, so that the underlying phase 
separation (PS) in the charge degrees of freedom should emerge at a macroscopic 
level, wiping out the SC properties. In fact, what we will see is that the SC pairs can 
survive even at non-zero magnetization, acquiring a modulation of FFLO type [7j. For 
higher values of the external field however, the density domain texture of the ground 
state emerges, and the SC pairs become first confined to the low density domains 
only, while the high density domains behave as polarized walls. Finally, above an 
upper value of the magnetic field superconductivity disappears and the ferromagnetic 
domains localize: the system becomes a ferrimagnetic insulator. 

In section II we describe some features of the model Hamiltonian, investigating 
in particular the role of the kinetic energy as to the onset of the different transitions. 
In section HI we investigate the static spin and charge structure factor dependence on 
the magnetization. The study of the Luttinger charge (Kp) and spin (i^cr) parameters 
allows to identify the presence of superconducting metallic and insulating phase. 
In section IV we explore the spatial dependence of density and magnetization with 
increasing the magnetic field. In section V we study the dependence of various type of 
pair-pair correlations and give evidence of the related emergence of FFLO oscillations. 
Finally in section VI we discuss the results and give some conclusions. 

2. The bond-charge Hubbard model 

2.1. Phase diagram at h = 

The ground- state phase diagram of Hamiltonian ([!]) at = 0, as obtained numerically 
in fT4| \TE\ [T6l U^9\ is reported in figure [TJ On the left-hand side this is given at fixed 
filling (n = 1) and varying u > 0, x, whereas on the right-hand side it is given at 
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Figure 1. Ground state phase diagram of H in one dimension (in the absence of 
external magnetic field). INS denotes the insulating phase, LL denotes the phase 
described by a Luttinger Liquid while BOW denotes the phase supporting bond 
order (charge density) wave [21 1 . Symbols (with joining dots) denote numerical 
results. Left panel: half-filling; continuous lines denote the analytical expressions 
for the critical lines obtained in [19 . Right panel: u = 1 (the legend refers to the 
methods used in [16] to locate the transition points). Dashed line marks x = 2/3. 

u = 1, varying n and x. Various phases can be recognized: in particular a Luther- 
Emery Uquid phase [21j with dominant superconducting correlations (SC) takes place 
in a range of filling values which depends on x. Within such phase, one can distinguish 
three different regimes: 

• (i)For < X < 1/2 the transition to the SC phase possibly takes place for n > 1 
and not too large u as soon as x 7^ 0, in agreement with bosonization predictions 
|22| (not shown). 

• (ii) For 1/2 < x < 2/3 we are in an intermediate regime. The SC phase still 
appears for < < Uc{x^ n) in a wider range ni < n < 2 with ni < 1. 

• (iii)For 2/3 < x < 1 nanoscale PS (NFS) is observed in the SC phase as a texture 
of two fluids of different densities ni and n^, which coexist up to Uc{x^n) for a 
range of filling values ni < n < rih < 2, with ni < 1. 

2.2. Role of kinetic energy 

The crossover among the above different regimes can be observed for instance on the 
ground-state energy e^^, plotted in figure [2f left) at half- filling and = versus x. 
With increasing x, egs turns out to increase linearly from the value cj^i^x — 0) — — 4/7r 
(non interacting electrons with spin) to the value Ck ~ — 2/7r reached for x = 1/2, 
in agreement with bosonization predictions. Being the latter coincident with the 
energy of a system of spinless fermions at quarter filling, from this point on the 
system becomes unstable with respect to PS into two regions at different densities 
(approximately 1/4 and 3/4) in region (ii), and correspondingly Cgs is seen to deviate 
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Figure 2. The case of vanishing magnetization at half-fiUing, u = 0. Left: 
Ground-state energy vs x. Right: the three contribution to kinetic energy in (|2]) 
plotted vs X. In both cases the three different regions (with white, gray, and dotted 
background) correspond to the three different behaviour described in section ^ 
Hubbard-hke, SC, NFS. 



from the linear increasing, reaching its maximum for x ~ 2/3. Above such value, the 
system enters region (iii), and the ground-state energy begins to decrease, to gain the 
value — 2/7r precisely at x = 1. At this point the system can be regarded either as 
describing a fluid of spinless fermions moving in a background of empty and doubly 
occupied sites [23], or as two coexisting fluids of single electrons moving in a background 
of empty and doubly occupied sites [17], in which case = 0.5 and rih = 1.5. 

The signiflcance of the points x = 1/2 and x = 2/3ati^ = becomes evident 
when representing the Hamiltonian in terms of on-site Hubbard projectors, which are 
operators deflned as = \a)i{P\i. Here \a)i are the states allowed at a given 

site i, and a = {0,t,|, 2} [|2) = | ti)]- this case the nonvanishing entries of the 
Hamiltonian matrix representation are read directly as the nonvanishing coeflicients 
of the projection operators. When rewritten in terms of these operators, H turns out 
to be a subcase of the more general Hamiltonian introduced by Simon and Aligia,[24j 

Hbc = - Y1 [^f + t.Xf-Xf + s^iXfX'f - Xf-Xf )] +u J2 

<ij>(T i 

in which tx = I — 2x and Sx = I — x. Besides the behaviour of Hbc is determined 
by the strength of tx and Sx. In particular, x = 1/2 implies tx = 0^ whereas below 
and above such value tx changes sign. A negative tx induces frustration in the motion 
of pairs, since it favors the presence in the ground state of momenta close to zero. A 
positive Sx term instead drives the pairs hopping favoring the modes with momenta 
close to the edge of the Brillouin zone. Hence we expect that for \tx\ ^ \sx\ -in our 
case X > 2/3- the mobility of the pairs becomes favored in the system for u < Uc{x). 
This is summarized in the right-hand side of flgure [2l where the three contributions 
to the kinetic energy in ([2]) are plotted separately at = 0. It is seen that in the 
three regions the relative weights of these contributions to the total kinetic energy are 
ordered in different ways. In particular, while at x < 1/2 the term with coefficient Sx, 
which does not conserve the number of pairs, plays the dominant role, for x > 1/2 this 
role is taken by the term describing the mobility of single electrons in a background 
of empty sites (coefficient 1), conserving the number of doublons; for x > 2/3 the 
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term with coefficient becomes the smallest, to vanish exactly at x = 1. In practice, 
a consequence of the above behaviour is that for x > 2/3 in the ground state the 
number of pairs is conserved, and the fluid behaves as a system of Ng single electrons 
moving in a background of L — Ns empty and doubly occupied sites. This is shown for 
instance in the charge structure factor N{q) = ^.((n^n^+^) — {ui) {uij^r))^^^^ ^ which 
correspondingly manifests a feature q = 27r Ng/L 

The reasons for the further choice of NPS, with respect to macroscopic PS, are 
to be found in the behaviour of the spin degrees of freedom, and their interplay 
with charge degrees of freedom. Indeed the term with coefficient Sx tends to favor 
antiferromagnetism at short range in an appropriate background of empty and doubly 
occupied sites. It is the surprising result of a very recent paper [19] that this simple 
assumption allows to map the Hamitonian in the SC regime into an effective XY model 
in a transverse field, recovering the critical line shown in figure [1] as the factorization 
line of the XY model. Remarkably, such line is not critical in the XY model: it is just 
the coupling of it to the charge degrees of in the original model which determines a 
change of phase. 

2.3. Results at h ^ 

At /i 7^ there are results available both at x = l[25], which case can be treated exactly 
also at non-zero temperature, and at x 7^ li26j, where the numerical method employed 
required the presence of a non-vanishing temperature. A low temperature peak in the 
specific heat, to be ascribed to the excitation of the spin degrees of freedom, is observed 
for X < 1/2 already at vanishing field, and at non- vanishing field for x < 1. 
Also, numerical results at T = have been obtained more recently for imbalanced 
species of ultracold fermionic atoms, both in case of attractive [27J, and in case of 
repulsive [18j ii. In this latter case it was realized that a non- trivial consequence of 
the phase coexistence in region (iii) is that within the NPS phase the number of pairs 
rid is constant with increasing m, up to a critical magnetization rUc = n — n^ = n — 2n(i 
at which the doublons begin to break, entering a regime of breached pairs. This is 
shown by the numerical data reported in figure [3l where also the analytical value 
obtained at the same tx values assuming Sx =^ are reported. It is seen that while in 
region (iii) the numerical data at = coincide with the theoretical curve already at 
m = 0, in the intermediate region (ii) this happens only at large enough values of m 
{m > rric). Also the size of the coexisting domains, microscopic for = 0, is observed 
to become macroscopic at appropriate non-zero magnetization in both regions. 

2.4- DMRG simulations 

Some observables considered in this work are better studied with periodic boundary 
conditions, while others are obtained numerically with more precision employing open 
boundary conditions. In both cases we retained up to 768 optimized DMRG states, 
performing usually three finite-system sweeps to enforce ground-state convergence. 
With more than 100 sites, the truncation errors in the density matrix weight remain 
O(10~^) or smaller. As far as the dependence on the system size is concerned, due 
to the large number of points selected in the phase diagram we decided to report the 
results for a fixed chain length L=160, unless otherwise specified. However, for some 
selected points we carried out a preliminar analysis of system size dependence (not 
reported here) and observed that the behaviour of the various correlation functions 




Figure 3. Number of doubly occupied sites versus magnetization at two different 
values of x, inside region (ii) (pink dots) and (iii) (violet dots) respectively. 
Continuous line is the theoretical prediction at x = 0.8 (sx = 0); dashed line 
is the critical magnetization at which the breached pair regime is entered in that 
case. 



was essentially the same at different values of (sufficiently large) L. Moreover, due 
to the onset of PS for some parameters value (see below), it could be difficult to 
reproduce the same qualitative spatial pattern by varying the system size. Hence the 
data presented at fixed length are chosen to be representative of the physical behaviour 
seen also at smaller sizes. 

3. Luttinger exponents: charge and spin structure factor at h ^ 

We have seen that at non-zero magnetization the structure of the ground state in 
region (iii) does not change -for what concerns the presence of empty and doubly 
occupied sites- up to the value rric. To some extent, this observation holds even inside 
region (ii). This is consistent with the possibility that an external magnetic field 
simply changes the polarization of the single electrons in the ground state. At non 
zero magnetization however the electrons should arrange differently regarding their 
pairing property. On general grounds, one of the two following scenarios is expected 
to hold above a first critical field at which the system begins to magnetize: either 
the SC properties are lost, or the SC pairs acquire a non- vanishing momentum and 
FFLO oscillations are observed [7J, due to the presence of polarized single electrons 
which do not pair. 

In order to explore what happens to the SC properties of the ground state 
in the two regions (ii) and (iii) at non-zero magnetization, we first evaluate the 
static charge N{q) and spin S{q) {S{q) = Ei((4^^4+l) - with 
s- = {ui^ — riii)) structure factors at different values of the average magnetization 
m = 1/(2L) < 5^^^ >. From the low frequency behaviour of these quantities one 
can extract the spin (K^) and charge (Kp) exponents which characterize the possible 
presence of gapped phases as: 

Kp = -N{q^O) , K, = -S{q^O) . 
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Figure 4. Luttinger exponent (see text) Kp u = 0, x = 0.6 (region (ii), left), 
and X = 0.8 (region (iii), right) as obtained from DMRG simulations (L=160). 
The two discontinuities identify the critical values mi and m2 discussed in the 
text. 

0.5 r 



0.4 




Figure 5. Spin structure factor at high magnetization m > rric (m = 0.4), 
for X = 0.4 and x = 0.6 as obtained from DMRG simulations (L = 160). The 
crossover from a spinful to a spinless regime is shown: continuous lines represent 
the theoretical prediction for a spinful fluid at half fllling (x = 0.4); and a spinless 
fluid of 2m fermions (x = 0.6). 



We recall that at m = both in region (ii) and (iii) the ground-state is a Luther-Emery 
SC liquid, i.e. it is characterized by K^j = and Kp > 1. The actual values of Kp for 
m > , as obtained from DMRG simulations at two different representative values of 
X, are given in figures 2] and [5] (the latter showing indirectly Kcr through the slope at 
^ = 0). 

At low enough magnetization, in both regions we still see evidence of a closed 
charge gap, and an open spin gap (not shown). However, as soon as m / in region 
(iii) Kp < 1 : the SC correlations are no longer dominant. On the contrary, in region 
(ii) Kp > 1 up to the critical magnetization mi ^0.1 reported in figure |H In this case 
we can conclude that superconductivity is still present in the magnetized ground state. 
Accordingly to our previous discussion, this should imply the presence of FFLO type 
of oscillations, which will be investigated in section [5l A further interesting aspect 
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Figure 6. Spin structure factor at various values of magnetization m < m2- Left 
panel: S(q) at x = 0.8, vertical lines in correspondence of ^Ugi- Right panel: 
S(q)/q at X = 0.6, vertical lines in correspondence of 2m . 



which appears in both regions is the presence of a higher magnetization value m2 above 
which the system appears to become insulating: even at finite size, Kp is zero within 
the numerical error. On the other hand our numerical data suggest that within the 
insulating regime the spin gap remains open until the magnetization reaches the value 
rric at which the breached pair regime is entered. Whereas for m > rric i^a = 1/2 
at variance with the standard Hubbard case in which K^j = 1 (see figure [5]). The 
result is a signal of the fully polarized nature of the single electrons for m > rric in 
regions (ii) and (iii), which implies that Kcr coincides with the Luttinger exponent of a 
gapless spinless liquid, i.e. 1/2. This is also confirmed by the fact that S{q) manifests 
a feature at g = Trng, where ris namely the number of single electrons, coincides for a 
fully polarize liquid with — n^. 

To exploit the nature of the difference between region (ii) and (iii) we also show in 
figure [6] 6'(g') for x = 0.6, 0.8 at various magnetization values below m2. We recall that 
in the standard Luttinger liquid case a feature is expected in S{q) for q — 27rncr[29j, 
whereas a signal of the presence of FFLO would be a feature at Q = 2i\(n^ — n|)[28]. 

Starting from region (iii), at x = 0.8 we observe a feature at = 27rn5|, which 
confirms that in this case the behaviour of magnetic correlations in the phase is that 
of a liquid of particle, moving in a background of n — empty and doubly 
occupied sites, in agreement with our previous observations: in this case, it is just 
the spin degrees of freedom of the single electrons which rearrange under the 
external magnetic field. In region (ii) instead, at x = 0.6 there is a neat feature 
at Q = 27r(n^ — n^), which is a further signal of the presence of FFLO oscillations. 
Such possibility will be explored in section [5l The feature is strongly reduced for 
m > mi. In both regions, the feature disappears at the value m2 discussed above. 

4. Phase separation and domain formation 

Due to quantum superposition and to the microscopic size of the domains, at vanishing 
magnetic field it is not possible to distinguish directly on the local density < n(j) > 
and magnetization < ^^(j) > profiles the presence of phase coexistence, even in region 
(iii) where it can easily be detected on the chemical potential [16j. The problem 
persists at low non-vanishing values of the magnetization, in which case however the 
local magnetization is observed to display a modulation with wavelength proportional 
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to either g or Q, depending on the value of x and consistently with the behaviour 
observed in S{q) (see previous section). With further increasing the magnetization, 
the presence of macroscopic phase separation (MPS) on both the local density and 
the local magnetization profiles appears well below the critical value rric at which the 
single electrons become fully polarized. This is seen in figures [3 and [H In fact, a 
more careful analysis shows that in region (ii) MPS appears as soon as m = mi (i.e. 
Kp = 1, where the system looses its SC properties) . Whereas in region (iii) MPS can 
be seen only above the magnetization m = m2 , in correspondence to the transition to 
the insulating state. As also shown from the profiles reported in figure [71 in this case 
the single electrons of the high density domains are fully polarized, and simultaneously 
the holes are expelled from the same domains, so that the system consists of localized 
domains with different average magnetization. 

In particular, in figure [3 we have chosen to plot separately in region (iii) the 
site dependence of the number of doubly < nd{j) > , singly occupied < nia > , 
and empty < n^(j) > sites in the different regimes. The presence of low and high- 
density domains can be easily recognized: these consists of spatial regions in which 
the single electrons move in a background of holes (low density) or doubly occupied 
(high density) sites, with different Fermi momenta (i.e., in this case different values of 
Us = — Usi^ shown in figure [7j). Since the single electrons can have different spin 
orientation only in the low density domains, the site magnetization turns out to be 
constant in the high-density ones, whereas it is modulated according to the number 
of single electrons with minority spins in the low density domains. The modulation 
disappears for m > mc, when the single electrons become fully polarized along all the 
system. A consequence of the spatial arrangement is that in the local magnetization 
profile the low and high density domains may change their relative spin orientation 
depending on the value of m. Indeed, as soon as m > m2 the system chooses to 
lower the magnetization of the low density domains, which weight is more relevant in 
the kinetic energy. Increasing m then amounts to align more and more spins of the 
electron in the low density domains, which at m = m2 are more numerous than those 
in the high density domains, so that the relative magnetization is reversed. Recalling 
that the high and low density domains are characterized respectively by the density 
of single electrons n^^, and n^/, in this case it is possible to calculate explicitly rric: 
rric = Usi - {rish - 1). 

As mentioned, in region (ii) it is possible to distinguish the presence of high and 
low density domains already at lower values of magnetization (figure [8]) , as soon as 
SC correlations cease to be dominant in the whole system (m = mi). In this case the 
single electrons with opposite spin orientation vary coherently with the sites in the 
low density domains, so that their difference (hence the local magnetization) remains 
constant there. We will see later that this is a further signal of the presence of FFLO 
in these domains. Whereas the local magnetization appears to be modulated in the 
high density domains, since the single electrons are not yet fully polarized. The low 
density phase is different from the one in region (iii), and the coherent behaviour 
of electrons with opposite spin suggest that in correspondence to the appearance of 
MPS SC pairs become confined to these domains. Further increasing m above m2 the 
behaviour resembles that of the x = 0.8 case: in both regions the system behaves de 
facto as an insulating ferrimagnet. 

Finally for m > rric the breached pairs regime is entered. Again, this fate is 
shared also by region (ii). 
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Figure 7. Site dependence of various densities for x = 0.8 (region (iii)) and 
1712 < m < rric ■ Left panel: local densities of holes and of doublons; right panel: 
local density of single electrons, and magnetization. 
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Figure 8. Site dependence of local densities of holes and doubly occupied sites 
(left panel); single electrons and magnetization (right panel) for mi < m < 1712 
ai X = 0.6. 



5. Pair-pair correlations and FFLO oscillations 

In this section we consider the possible presence of FFLO oscihation in pair pair 
correlations at non-zero magnetization via the study of various pairing correlators 
(AlAi-^R) where the operator Ai at site may refer to A^^^^^^ = — Q^Qt (o^site 
pairing) or to singlet and triplet combinations on adjacent sites i and z + 1 



A^l""^ = Si = — (ci+i;Ci^ - Ci^i^Cii) 

^trip _ _ _^ Ci-^i^Cii) 



^trip- 



(3) 



For each one of these operators, as well as for their correlators, we can use Hubbard 
operators X and give a decomposition based on the high-density subspace made only of 




Figure 9. Various types of correlators (A] with the choices for Ai reported 

in the legend (compare with (O). Considering the emerging PS we have selected 
two possible starting sites z = 45 or z = 60. In the left hand side the correlation 
functions are plotted in reciprocal space and the arrow marks the expected FFLO 
peak at gpFLO = ^Ft ~ ^Fl = 7T(n^ — n^) = 27Tmz. Here the magnetization 
density, that allows for the FFLO effect, is ruz = 0.05. The oscillating pattern in 
real space is evident from the right panel, especially for the dominating ?7-pairs 
correlation that decay slowly. The Ti correlations instead do not exhibit the same 
effect. The model parameters here are u = and x = 0.6. 



doubly occupied sites and singly occupied sites {2a elements), the low-density subspace 
where no doubly occupied sites exist while empty sites are permitted (Oa elements). 
Mixed contributions also must be considered, as can be seen from the example of 

T mixed vX-'^ v^T i v^T vX-'^ 

For T-ii it is sufficient to reverse the index a of single occupation. The component 
of the triplet and the singlet have similar, although longer, expressions that we do not 
report here for the sake of brevity. 

We discuss the results of a quantitative DMRG analysis by choosing x = 0.6 and 
X = 0.8 as representative parameters of the two somehow different situations of region 
(ii) and (iii) respectively. All the numerical data presented in this section are taken 
on chains of L = 120 with open boundary conditions (OBC) sites at half filling. First, 
in figure [9] we show evidence of FFLO behaviour showing up when the populations 
are slightly imbalanced in region (ii). In figure [TOl instead, we can see how the FFLO 
effect is progressively wiped out if either the magnetization density is too high (central 
panel), or if the interaction u is increased (right panel). 

Due to the possible emergence of PS, that is spatial inhomogeneities, it will be 
useful to inspect also the variation of the correlators {A^-Aij^R) along the chain at fixed 
as a function of i. When R = 1 there is a site i + 1 in common between and 
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Figure 10. Left and central panel: as figure [9] but with magnetization density 
increased to 0.1 and 0.15 respectively. In the former case a weak signature of 
FFLO peak remains (it would be expected in correspondence of the dotted vertical 
line), while in the latter case is wiped out. The same happens by increasing the 
Coulomb interaction at the same value of =0.1 (for instance w = 1 in the 
example of the right panel). 



AiJ^i (apart from the case Ai = rji) and there are some simphfications in the fohowing 
correlation functions referred to the "high" and "low" sectors of the Hilbert space 

where ± in the subscript corresponds to a =t, i and a means spin flip. In addition, 
it turns out that the "high" parts of the singlet and of the 0-component triplet are 
always equal except for the sign 

(:r^:'^^T^:t\} = -is't^'^'s^if) (4) 

In figure [11] we examine the local dependence of pairing correlations along a chain 
showing how PS appears also in SC properties. The relative distance is now fixed to 
R = 1 site and the leftmost site i is varied; again we select the example u = 1, x = 0.6 
as in figure [TOl (right). As long as the net magnetization density is small the dominant 
pairing correlations are of singlet type and the dependence on i is not strong. When 
the magnetization is increased the +1 component of the triplet starts to dominate, but 
in the high-density regions the pairing correlations are generally suppressed, with the 
exception of a small enhancement of the —1 component of the triplet. Notice that we 
may arrive at a situation in which there are no unpaired down particles (bottom right 
panel) and, correspondingly, the singlet and Tq correlators acquire the same profile 
except for the sign as in (j4]). 

Having identified the dominant type of pairing correlations, at least at short 
distance, we can also inspect the separate contributions from low- and high-density 
sectors of the Hilbert space to the singlet and T+i correlators, according to expressions 
above. This is done in figure [T2l 

Qualitatively the result can be summarized as follows. On increasing the 
magnetization m the spin-up electrons align along the field first in the high-density 
islands and then in the low-density ones. When in region (ii) the system exhibits 
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Figure 11. Various types of correlators (A^Ai^fi) with the choices for Ai written 
in ^ and the distance fixed to R = 1. The model parameters here are u = 1, 
X = 0.6 and with OBC one can appreciate the spatial dependence on i. Upper 
panel: ruz =0.1 (left) and rriz = 0.15 (right); lower panel: niz = 0.25 (left) and 
rriz = 0.40 (right). Due to refiection symmetry only the left half of the chain is 
displayed. 



FFLO oscillations in pair structure factor. When the domains finally appear at a 
macroscopic level (m = mi), the FFLO pairs become mostly localized into the low 
density domains, and the corresponding peak is weakened; globally the system behaves 
as metal made of ferromagnetic domains alternating with superconducting ones. The 
SC properties of the low density domains fully disappear only when single electron 
localize (m > m2), as shown by the right panel of figure [TTl 
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Figure 12. Same as figure [TT] for the singlet and T+i correlation functions at 
distance R = 1, where full up triangles now denote the "high" component and the 
empty down triangles denote the "low" one. Left panel is for =0.1 while the 
right one is for = 0.25. 

6. Discussion and Conclusions 

The results of the previous sections provide a coherent scenario of the physics of the 
SC phase described by the Hamiltonian ([1]) for x > 1/2 in presence of non- vanishing 
magnetization. With respect to the case x < 1/2 (region (i)) , the different behaviour 
observed previously at m = persists up to high magnetization values (m > mc), as a 
consequence of the rearrangement in the kinetic energies of holes and doublons which 
takes place at x = 0.5 (see section II). The main features of the phase diagram in the 
X — m plane are schematically shown in figure [13] for regions (ii) and (iii). 

At low magnetization the diagram shows evidence of the two regions, which turn 
out to behave differently up to m = m2; here m2 is the value at which the holes and 
minority spin single electrons localize into the low density domains, and the doublons 
localize into the high density domains: the system becomes insulating. In region (ii) 
SC persists in the whole system at non-zero magnetization up to m = mi < m2; the 
pair structure factor correspondingly shows a strong FFLO peak. Macroscopic phase 
separation appears on both the density and the magnetization profiles for m > mi: in 
this case SC correlations become mainly confined to the low-density domains, whereas 
the system is still metallic due to the mobility of the single electrons. The low and high 
density regions display different magnetizations, so that for mi < m < m2 the system 
can be considered an itinerant ferromagnet, with superconducting FFLO domains. As 
for region (iii), there is no evidence for dominant SC correlations in the NFS phase 
as soon as m 7^ 0; the system behaves as a liquid of single particles moving in a 
background of empty and doubly occupied sites textured into domains of microscopic 
size[16j. Macroscopic phase separation into high- and low-density domains appears 
on the density and magnetization profiles only for m > m2, simultaneously with the 
transition to the insulating state. Finally, at a higher magnetization mc, the regime 
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Figure 13. Schematic phase diagram at < w < Uc(x) and half-fiUing, in the 
X — m plane for x >l/2 

of breached pairs described in [18j is entered. 

We emphasize that the presence of a FFLO regime in the repulsive Hubbard model 
takes place for a range of values of the off-diagonal Coulomb repulsion (l/2<x<2/3) 
in principle accessible to experiments. The results reported here for such regime are 
complementary to those obtained in [30j for the extended Hubbard model; in fact, they 
are expected to persist also in presence of non-vanishing nearest neighbor Coulomb 
repulsion, since similar incommensurate spin and charge correlations are observed also 
in that case [31j. The results are reminiscent for some aspects of those reported in [32] 
for the standard Hubbard model in the attractive u region. We guess that our FFLO 
regime is to be identified to the weak LO (Larkin-Ovchinnikov) regime described there, 
whereas the regime of macroscopic phase separation of the high density polarized walls 
and the low density FFLO domains could well coincide with the strong LO regime 
discussed there, in which pair-pair correlations exhibit spatial dependence. The main 
difference we see with [32] is that in the present case with increasing the magnetization 
the domain walls (to be identified with our high density domains) tend to localize, 
in contrast to the delocalized behaviour observed there. Ultimately this feature is 
responsible for the transition to the insulating regime found here. We expect that the 
two scenarios could merge in the region of weakly attractive u and x < 0.5. 

As mentioned, with varying magnetization a second transition then occurs in both 
regions at m = m2; in this case the system consists of two segregated metals of fully 
polarized single electrons with different Fermi momenta, which as a whole behave as 
a (ferrimagnetic) insulator. In two dimensions, we argue that the behaviour of such 
insulator could well be that of metallic stripes; whereas for mi < m < m2 in region 
(ii) the physics chould be that of a striped superconductor [33j. 

The analysis is performed here at finite size. Due to the different qualitative 
behaviour of the phases depicted in figure [13] the scenario should be considered quite 
plausible. A finite size scaling analysis is expected to provide a quantitative accurate 
derivation of the critical lines. In particular, since our study just considered two 
far apart typical points {x = 0.6 for region (ii), and x = 0.8 for region (iii)), it is 
not possible to infer from it whether the change from region (ii) to (iii) is a smooth 
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crossover or a true transition line. 
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